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The MiniBooNE experiment at Fermilab reports results from an analysis of the combined Ve and 
Ve appearance data from 6.46 x lO^'' protons on target in neutrino mode and 11.27 x 10^° protons 
on target in antineutrino mode. A total excess of 240.3 ± 34.5 ± 52.6 events (3.8o") is observed from 
combining the two data sets in the energy range 200 < < 1250 MeV. In a combined fit for 

CP-conserving Ve and v^^ Ve. oscillations via a two-neutrino model, the background-only 

fit has a x^-probability of 0.03% relative to the best oscillation fit. The data are consistent with 
neutrino oscillations in the 0.01 < Am^ < 1.0 eV^ range and with the evidence for antineutrino 
oscillations from the Liquid Scintillator Neutrino Detector (LSND). 



There is growing evidence for short-baseline neutrino 
anomalies occuring at an L/En ~ 1 m/MeV, where is 
the neutrino energy and L is the distance that the neu- 
trino travelled before detection. These anomalies include 
the excess of events observed by the LSND [l| and Mini- 
BooNE 0-13 experiments and the deficit of events ob- 
served by reactor Q and radioactive-source experiments 
0. There have been several attempts to interpret these 
anomalies in terms of 3-l-N neutrino oscillation models 
involving three active neutrinos and N additional ster- 
ile neutrinos [7-12[. (Other more exotic explanations 
include, for example, Lorentz violation [isj and sterile 
neutrino decay [1J|.) This paper presents a combined 
oscillation analysis of the MiniBooNE Ve. and appear- 
ance data, corresponding to 6.46 x 10^" protons on tar- 
get (POT) in neutrino mode and 11.27 x 10^" POT 
in antineutrino mode, which is approximately twice the 
antineutrino data reported previously 



This analysis fits both z^^ 



and 



oscilla- 



tions with the same oscillation model over the full neu- 
trino energy range 200 < £'^^ < 3000 MeV, where E^^ 



is the reconstructed neutrino energy assuming quasielas- 
tic scattering kinematics |15| . The neutrino oscillation 
energy region is defined to be 200 < £'^^ < 1250 MeV, 
which is where an LSND-like signal (same L/E^) is ex- 
pected. Combining neutrino and antineutrino data over 
the full energy range has the advantage of decreasing sta- 
tistical and systematic errors. The analysis assumes no 



significant 



or Vf. disappearance. This simpli- 



fication may change the fitted — > v,. and z/^ — ^ Ve 
appearance probabilities by up to ~ 20%. Furthermore, 
it has been suggested that nuclear effects associated with 
neutrino interactions on carbon can affect the reconstruc- 
tion of the neutrino energy and the determination of the 
neutrino oscillation parameters [l6| . These effects are 
not fully accounted for in the analysis and may affect 
somewhat the oscillation fit parameters discussed below. 

The neutrino (antineutrino) flux is produced by 8 GeV 
protons from the Fermilab Booster interacting on a beryl- 
lium target inside a magnetic focusing horn set at posi- 
tive (negative) polarity. In neutrino (antineutrino) mode, 
positively (negatively) charged mesons produced in p-Be 



2 



interactions are focused in the forward direction and sub- 
sequently decay primarily into (i'^)- The flux of neu- 
trinos and antineutrinos of all flavors is simulated using 
information from external measurements [l7j . In neu- 
trino mode, the v^, D^, v^, and Dg flux contributions at 
the detector are 93.5%, 5.9%, 0.5%, and 0.1%, respec- 
tively. In antineutrino mode, the Vf^, Dg, and h'g flux 
contributions at the detector are 83.7%, 15.7%, 0.4%, 
and 0.2%, respectively. The t'^ and i?^ fluxes peak at 
approximately 600 MeV and 400 MeV, respectively. 

The MiniBooNE detector is described in detail in ref- 
erence [S]. The detector is located 541 m from the 
beryllium target and consists of a 40- foot diameter sphere 
filled with 806 tons of pure mineral oil (CII2). Neutrino 
interactions in the detector produce charged particles 
(electrons, muons, protons, pions, and kaons) which in 
turn produce scintillation and Cherenkov light detected 
by the 1520 8-inch photomultiplier tubes (PMTs) that 
line the interior of the detector and an optically iso- 
lated outer veto region. Event reconstruction and par- 
ticle identification are derived from the hit PMT charge 
and time information. 



The signature of — ^ and f^, 



oscillations 



is an excess of Ve and i^g-induced charged-current quasi- 
elastic (CCQE) events. Reconstruction [l^ and selection 
requirements of these events are almost identical to those 
from previous analyses 0, 0] with an average reconstruc- 
tion efficiency of ^ 10 — 15% for events generated over the 
entire volume of the detector. Recent improvements to 
the analysis include a better determination of the intrin- 
sic Ve background from K'^ decay through the measure- 
ment of high-energy neutrino events in the SciBooNE 
experiment (20l |. a combined error matrix for neutrino 
and antineutrino data with correlated and uncorrelatcd 
errors, a better determination of neutral-current tt^ and 
external event background in antineutrino mode due to 
the increase in statistics of the antineutrino mode data 
sample, and the use of a likelihood fit with frequentist 
corrections from fake data studies for both the neutrino- 
mode and antineutrino-mode data. The detector cannot 
distinguish between neutrino and antineutrino interac- 
tions on an cvent-by-event basis. However, the fraction 
of CCQE events in antineutrino (neutrino) mode that are 
due to wrong-sign neutrino (antineutrino) events was de- 
termined from the angular distributions of muons created 
in CCQE interactions and by measuring charged-current 
single TT^ 

and De CCQE background events 



events [21 
The predicted Ve 
for the neutrino oscillation energy range 200 < E^^ < 
1250 MeV are shown in Table [J for both neutrino mode 
and antineutrino mode. The predicted backgrounds to 
the Ve and Dg CCQE sample are constrained by measure- 
ments in MiniBooNE and include neutral current (NC) 
7r° events [l^l with photonuclear interactions, A — > iV7 
radiative decays psj . and neutrino interactions external 
to the detector. Other backgrounds from mis-idcntificd 



TABLE I: The expected (unconstrained) number of events for 
the 200 < -E?^ < 1250 MeV neutrino oscillation energy range 
from all of the backgrounds in the and appearance anal- 
ysis and for the LSND expectation of 0.26% oscillation prob- 
ability averaged over neutrino energy for both neutrino mode 
and antineutrino mode. 



Process 


Neutrino Mode Antineutrino Mode 


v,, & u„ CCOE 


37.1 


12.9 


NC 7r° 


252.3 


112.3 


NC A ^ iV7 


86.8 


34.7 


External Events 


35.3 


15.3 


Other 1^,1 & 


45.1 


22.3 


Ve & Ve from /i* Decay 


214.0 


91.4 


Ve & Ve from Decay 


96.7 


51.2 


Ve & Ve from Decay 


27.4 


51.4 


Other Ve & Ve 


3.0 


6.7 


Total Background 


797.7 


398.2 


0.26% V^ Ve 


233.0 


100.0 



or v^ [2J, |25[ and from intrinsic Vg and Vg events from 



the TT ^ n decay chain arc constrained and obtain their 
normalizations from the Vf^ and Dp, CCQE data samples, 
which consist of 115,467 (50,456) events in neutrino (an- 
tineutrino) mode in the 200 < Ef^^ < 1900 MeV energy 
range. 

Systematic uncertainties are determined by consider- 
ing the predicted effects on the i'^, P^, Vg, and Vg CCQE 
rate from variations of actual parameters. These include 
uncertainties in the neutrino and antineutrino flux esti- 
mates, uncertainties in neutrino cross sections, most of 
which are determined by in situ cross-section measure- 
ments at MiniBooNE, and uncertainties in detector mod- 
eling and reconstruction. A covariance matrix in bins 
of E^^ is constructed by considering the variation from 
each source of systematic uncertainty on the Vg and Vg 
CCQE signal, background, and i/^ and CCQE predic- 
tion as a function of E^^ . This matrix includes corre- 
lations between any of the Vg and Dg CCQE signal and 
background and Vfj_ and CCQE samples, and is used 
in the calculation of the oscillation fit. 

Fig. [T]shows the E^^ distribution for Vg and i?g CCQE 
data and background in neutrino and antineutrino mode 
over the full available energy range. Each bin of re- 
constructed E^^ corresponds to a distribution of "true" 
generated neutrino energies, which can overlap adjacent 

TABLE II: The number of data, fitted (constrained) back- 
ground, and excess events in the Ve and Ve analyses for neu- 
trino mode, antineutrino mode, and combined in the neutrino 
oscillation energy range 200 < E^^ < 1250 MeV. The un- 
certainties include both statistical and constrained systematic 
errors. All known systematic errors are included in the sys- 
tematic error estimate. 



Mode 


Data 


Background 


Excess 


Neutrino Mode 


952 


790.0 ±28.1 ±38.7 


162.0 ±47.8 


Antineutrino Mode 


478 


399.6 ± 20.0 ± 20.3 


78.4 ± 28.5 


Combined 


1430 


1189.7 ±34.5 ±52.6 


240.3 ± 62.9 



3 



Neutrino 

• Data (Stat err.) 
I I Ve from fi*'" 
I I Ve from K*'" 
I I V from K° 

X misid 
I I A ^ My 
H| dirt 

other 

Constr. Syst. Error 




0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
0.4 
0.3 
0.2 
0.1 
0.0 



-0.1 



0.2 



Data - expected background 
sin^29=0.004^m2=1 .OeV^ 
sin^29=0.2,Am2=0.1eV^ 
Combined 2v Best Fit 
Combined 3+2 Best Fit 

Neutrino 




Antineutrino 



0.4 
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E?^/GeV 



FIG. 1: The neutrino mode (top) and antineutrino mode (bot- 
tom) E^^ distributions for Ue CCQE data (points with sta- 
tistical errors) and background (histogram with systematic 
errors) . 

bins. In neutrino (antineutrino) mode, a total of 952 
(478) events pass the Ve event selection requirements with 
200 < E^^ < 1250 MeV, compared to an expectation of 
790.0 ±28.1 ±38.7 (399.6 ±20.0 ±20.3) events, where the 
first error is statistical and the second error is systematic. 
This corresponds to a neutrino (antineutrino) excess of 
162.0 ± 47.8 (78.4 ± 28.5) events. Combining the data in 
neutrino mode and antineutrino mode, the total excess 
is 240.3 ±62.9 events. Fig. [2] shows the event excesses 
as a function of E^^ in both neutrino and antineutrino 
modes. The number of data, fitted background, and ex- 
cess events for neutrino mode, antineutrino mode, and 
combined arc summarized in Table HIl 

Many checks have been performed on the data, includ- 
ing beam and detector stability checks that show that 
the neutrino event rates arc stable to < 2% and that 
the detector energy response is stable to < 1% over the 
entire run. In addition, the fractions of neutrino and an- 
tineutrino events are stable over energy and time, and 
the inferred external event rate corrections are similar in 
both neutrino and antineutrino modes. 

A comparison between the MiniBooNE and LSND an- 
tineutrino data sets is given in Fig. [3l which shows the 
oscillation probability as a function of L/ E^ for v^^ — > i^e 
and — >■ Df. candidate events in the L/ E^ range where 



FIG. 2: The neutrino mode (top) and antineutrino mode (bot- 
tom) event excesses as a function of E^^ . Also shown are the 
expectations from the best two-neutrino and 3-1-2 joint oscilla- 
tion fits with 200 < S?-^ < 3000 MeV and from two reference 
values in the LSND allowed region. All known systematic er- 
rors are included in the systematic error estimate. 

MiniBooNE and LSND overlap. The data used for LSND 
and MiniBooNE correspond to 20 < E^ < 60 MeV and 
200 < E^^ < 3000 MeV, respectively. The oscilla- 
tion probability is defined as the event excess divided 
by the number of events expected for 100% Vg 
and — > De transmutation in each bin, while L is the 
distance travelled by the neutrino or antineutrino from 
the mean neutrino production point to the detector and 
Ejj is the reconstructed neutrino or antineutrino energy. 
The largest oscillation probabilities from both LSND and 
MiniBooNE occur at L/E^ > 1 m/MeV. 

The MiniBooNE data are next fit to a two-neutrino 
oscillation model, where the probability, P, of i^,, — ^ 
h'e and i>i^ — > Ve oscillations is given by P = 
81X1^ 298111^ {1.27 Am^L/E^), sin^ 26^ = 4|[/e4p|C/,.4p, and 
Am^ = Amli = rn\ — m\. The oscillation parameters 
are extracted from a combined fit to the Ve, v^, v^, and 
Dfj, CCQE event distributions. The fit assumes CP con- 
servation with the same oscillation probability for neu- 
trinos and antineutrinos, including both right-sign and 
wrong-sign neutrinos, and no significant z/^, D^, Ve^ov 
disappearance. Using a likelihood-ratio technique [J], the 
best oscillation fit for 200 < S,?^ < 3000 MeV occurs at 
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FIG. 3: The oscillation probability as a function of L/E^^ 
for and — > i/g candidate events from MiniBooNE 

and candidate events from LSND. The data points 

include both statistical and systematic errors. Also shown are 
the oscillation probabilities from the two-neutrino and 3-f 2 
joint oscillation fits. 

(Am^, sin^26l) = (0.037 cV^, 1.00). The ^ Indf for the 
best-fit point in the neutrino oscillation energy range of 
200 < E^^ < 1250 MeV is 24.7/15.6, corresponding to a 
probability of 6.7%. The probability of the background- 
only fit relative to the best oscillation fit is 0.03%. Fig.|3] 
shows the MiniBooNE closed contours for and ap- 
pearance oscillations in neutrino mode and antineutrino 
mode separately in the 200 < E^^ < 3000 MeV en- 
ergy range, where a two-neutrino oscillation model is 
assumed and where frequentist studies were performed 
to determine the confidence level (C.L.) regions. The 
separate best fits for neutrino mode and antineutrino 
mode arc at (Am^, sin^ 26*) values of (3.14 eV^, 0.002) 
and (0.05 eV^, 0.842). In the neutrino oscillation en- 
ergy range of 200 < £'^^ < 1250 MeV, the x'^/ndf 
for the best-fit points in neutrino mode and antineu- 
trino mode are 13.2/6.8 and 4.8/6.9 with probabilities 
of 6.1% and 67.5%, respectively. The background-only 
fit has a x^-probability of 1.6% and 0.5% relative to the 
best oscillation fits for neutrino and antineutrino, respec- 
tively. Fig. [5] shows the closed contours for the com- 
bined fit. The allowed regions for Ato^ < 1 eV^ are in 
agreement with the LSND allowed region [J and con- 
sistent with the limits from the KARMEN experiment 
[iil l . Fig. [2] shows the expectations from both the best 
two-neutrino joint oscillation fit and from a 3-1-2 joint os- 
cillation fit as a function of E^^ in both neutrino and 
antineutrino modes. The best-fit parameters from the 
3-^2 oscillation fit are Amj-^ = 0.082 eV^, Am^^ = 0.476 
eV2, |C/e4p|C/^4|' = 0.1844, pesl^p^^? = 0.00547, and 
(j) = I.OOOStt. The 3-1-2 fit has three more parameters 
than the two-neutrino fit 12| and will be discussed in a 
future publication. 

In summary, the MiniBooNE experiment observes a 
total excess of 240.3 ± 62.9 I'e and events (3.8cr) in 
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FIG. 4: MiniBooNE allowed regions in neutrino mode (top) 
and antineutrino mode (bottom) for events with E^^ > 200 
MeV within a two-neutrino — > and Vii, — >■ oscillation 
model. Also shown is the i/^ — > limit from the KARMEN 
experiment [2^. The shaded areas show the 90% and 99% 
C.L. LSND i?^ Ue allowed regions. The black stars show 
the best fit points, while the crosses show LSND reference 
values. 
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FIG. 5: MiniBooNE allowed regions in combined neutrino and 
antineutrino mode for events with 200 < E^'^ < 3000 MeV 
within a two-neutrino i/^ — >■ i/e and i/^ — >■ i/e oscillation model. 
Also shown is the limit from the KARMEN exper- 

iment [H. The shaded areas show the 90% and 99% CL. 
LSND Vf_i — 7- allowed regions. The black star shows the 
best fit point. 

the neutrino oscillation energy range 200 < E^^ < 
1250 MeV. The allowed regions from a two-neutrino fit 
to the data, shown in Fig. [5l arc consistent with — > 
and — )■ Ve oscillations in the 0.01 to 1 eV^ Am^ range 
and consistent with the allowed region reported by the 
LSND experiment [ij. 
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